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ABSTRACT 

The  communication  channel  between  a  satellite  and  a  ship  earth  station  (SES)  is 
described  by  a  model  which  includes  multipath  fading,  doppler  shift  and  noise.  Multi- 
path  fading  is  caused  by  reflections  from  the  sea  surface.  These  reflections  can  affect  the 
system  performance,  especially  at  low  elevation  angles  or  when  SES  is  using  low  gain 
antennas.  Doppler  shift  is  a  ver\'  important  effect  when  using  low  altitude  satellites,  be- 
cause of  the  high  velocities  involved. 

This  thesis  describes  and  presents  a  software  simulator  for  multipath  fading  in  the 
maritime  communications  environment.  Analysis  of  throughput  of  an  unslotted  Aloha 
maritime  mobile  satellite  communication  channel  is  also  presented. 
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I.     INTRODUCTION 

Future  satellite  systems  for  maritime,  aeronautical  and  land  mobile  applications 
(with  mobile  terminals  for  low  data  rate  and  digital  voice  transmission)  will  have  small 
antennas  having  gains  between  2  and  12  dB  and  a  broad  beamwidth.  These  antennas 
pick,  up  multipath  reflections  of  satellite  signals  causing  severe  fading  of  the  received 
signals. 

Several  fading  simulators,  using  hardware,  software,  or  both,  have  been  designed. 
Among  them  are:  JPL,  CRC  and  Virginia  Tech.  All  of  them  deal  with  land  mobile 
communications  and  with  geostationary'  satellites. 

In  the  future,  low  altitude  satellites  will  become  important  ahernatives  for  mobile 
communications  because  of  the  crowded  geostationary'  orbit.  Studies  have  demonstrated 
that  five  satellites  in  circular  orbit  will  provide  continuous  single  coverage  of  the  earth; 
seven  satellites  will  provide  double  coverage  at  some  minimum  elevation  angle.  By  the 
same  criterion  it  was  found  that  eleven  satellites  can  provide  continuous  triple  coverage 
and  fourteen  satellites,  quadruple  coverage.  Continuous  single  coverage  means  that  any 
two  stations  on  earth  can  communicate  using  this  constellation  of  five  satellites  at  any 
moment.  Double  coverage,  on  the  other  hand,  will  provide  an  additional  communication 
path;  triple  coverage  will  provide  three  dilTerent  paths  and  so  on. 

Besides  the  economical  achievement  of  maximum  coverage,  many  other  factors  have 
to  be  taken  into  account  in  system  design.  These  results  may  prove  of  some  interest  to 
planners  of  satellite  systems  for  mobile  communications  and  surveillance  [Ref  1] 

Low  earth  orbits  (LEO)  are  generally  circular  or  slightly  eccentric.  In  order  to 
simphfy  the  analysis  and  equations  described  in  this  thesis,  the  focus  will  be  on  circular 
orbits. 

This  thesis  presents  a  fading  simulation  program  for  mobile  maritime  satellite  com- 
munications involving  low  altitude  satellites. 

Chapter  II  describes  the  different  signals  involved  in  the  maritime  communications 
channel,  the  atmospheric  effects  that  affect  them,  and  also  demonstrates  the  application 
of  the  theor}"  of  Rician  fading  statistics. 

Chapter  III  describes  the  simulation  program,  explains  the  purpose  of  each  of  its 
components  and  illustrates  some  results  obtained  by  running  several  simulations  for 
different  values  of  frequency  and  satellite's  velocity. 


Chapter  IV  describes  the  throughput  of  unslotted  Aloha  with  power  capture  prob- 
ability in  a  Rician  fading  environment. 

Finally,  Appendix  A  contains  a  listing  of  the  computer  program  of  the  simulator, 
and  Appendix  B  shows  the  derivation  of  the  formulas  used  in  chapter  IV. 


II.     PHYSICS  OF  MARITIME  MOBILE  SATELLITE  PROPAGATION 


A.  INTRODUCTION 

Any  signal  in  maritime  mobile  communications  is  typically  dominated  by  a  multi- 
path  fading  process  caused  by  reflections  coming  from  the  surface  of  the  sea.  These 
multipath  signals  are  expected  to  be  stronger  at  low  satellite  elevation  angles,  where  the 
main  lobe  of  the  antenna  of  the  ship  earth  station  (SES)  picks  them  up.  No  specular 
component  is  expected,  except  for  the  very  calm,  mirror-like  sea;  therefore,  a  model  with 
Rician  statistics  is  suggested  in  this  chapter. 

B.  MULTIPATH  FADING  AMPLITUDE 

It  is  well  known  that  multipath  rellections  from  the  sea  surface  consist  of  specular 
and  diffuse  components.  The  relative  strength  of  these  components  is  due  to  wave  height 
and  to  elevation  angle,  with  the  former  the  dominant  factor  [Ref  2].  For  a  smooth  sea 
state  the  specular  component  dominates.  As  wave  height  increases,  the  diffuse  compo- 
nent also  increases,  and  over  a  wave  height  of  0.5  m,  the  diffuse  component  becomes 
constant.  Sea  state  probability  studies  have  been  done  in  the  past  and  it  was  concluded 
that  wave  heights  are  usually  in  a  0.5  to  5  m.  range  [Ref  3].  Based  on  these  observations, 
it  is  assumed  that  sea  conditions  will  be  rough  most  of  the  time  and  thus  amplitude 
fading  can  be  expressed  in  statistical  terms. 

C.  STATISTICAL  MODEL  FOR  MULTIPATH  FADING  AMPLITUDE 

For  a  satellite  maritime  link,  the  dominant  multipath  component  is  the  diffuse 
component.  The  signal  can  therefore  be  represented  by  two  components:  direct  and  dif- 
fuse (see  Figure  1). 


DIRECT  WAVt 


Figure  1.      Propagation  model  for  maritime  mobile  satellite  commmiications 

The  direct  component  is  the  line  of  sight  (LOS)  signal  from  the  satellite.  The  dilTuse 
component  is  the  resultant  of  many  random  scattered  waves  arriving  in  a  uniform  an- 
gular distribution  [Ref  4].  If  each  of  the  waves  is  represented  by  a  random  phasor,  the 
total  diffuse  component  at  the  output  of  the  antenna  is 


(2T) 


(=1 


where 

•  r  is  the  amplitude  of  the  difliise  component 

•  0  is  the  phase  of  the  diffuse  component  with  respect  to  the  direct  component 


•  A,  is  the  amplitude  of  the  i"'  scattered  wave  with  respect  to  the  direct  component. 

If  the  scattered  waves  are  sufficiently  random  and  their  phase  is  uniformly  distrib- 
uted over  2n,  Beckmann  shows  that  the  diffuse  component  will  be  Rayleigh  distributed 
with  a   density  function  given  by  [Ref  5] 

/.W  =  ^exp[-4]  ''^0  (2.2) 

p{r)  =  0  r  <  0 

where 

•  r  is  the  received  diffuse  signal  level  (voltage) 

•  a  =  £{H} 

The  Rayleigh  distribution  requires  the  scattered  waves  to  have  a  uniform  phase  dis- 
tribution over  an  interval  of  2;r.  If  the  terrain  does  not  provide  for  this  requirement,  the 
Rayleigh  distribution  will  not  be  valid  [Ref  6]. 

The  vector  sum  of  a  direct  signal,  with  a  Rayleigh  fading  component,  results  in  a 
composite  signal  with  Rician  envelope  statistics  [Ref  7],  and  the  probability  density 
function  is: 

/.(/-)  =  i?-exp    -— — ^    /o(-r^)         r>0  (2.3) 


p{r)  =  0  r  <  0 

where 

•  r  represents  the  envelope  of  the  Rayleigh  fading  signal 

•  a  =  E{r'] 

•  Aq  is  the  amplitude  of  the  direct  signal 

•  /q  is  the  modified  Bessel  function  of  zero  order 

By  normalizing  the  direct  component  to  unity,  this  reduces  to: 

piR)  =  M.,J_ll±±y^^2K^      R^o  (2.4) 


P{R)  =  0 


R<0 


where 

•      A  is  equal  to  the  power  in  the  multipath  component  over  the  power  in  the  direct 
component,  that  is: 

a 


k  = 


Q 


•      R  is  the  ratio  of  the  amplitude  of  the  fading  signal  to  the  amplitude  of  the  direct 
component 

The  parameter  k  is  called  the  Rice  factor.  It  is  considered  a  key  parameter  in  the 
Rician  distribution,  and  its  decibel  equivalent  is  K  =  10  log,oA-.  Figure  2  illustrates  se- 
veral Rician  distributions  for  different  values  of  A„. 


The  phase  of  the  received  signal  is  not  uniformly  distributed  [Ref  5]  and  its  density 
function  is  given  by: 
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p{e)  =  -^txpi--f-l]  +  G^^ntxp{G'){l+erAG))]      0<e<2n         (2.5) 


where 


•    G  is  defined  as: 
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•     erf  (G)  is  the  error  function  defmed  as: 
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Figure  2.      Rician  density  function  for  a  =   1,  f  =  869  MHz,  V  =  0  Kts. 

1.     Cumulative  Distribution  Function  (CDF) 

The  cumulative  distribution  function  (CDF)  refers  either  to  the  function  G(A) 
for  which  the  signal  exceeds  a  level  A,  or  to  the  function  F(A)  for  which  the  signal  is 
below  a  level  A.  The  former  is  the  one  used  in  this  thesis. 

The  theoretical  expression  for  the  Rician  CDF  is  obtained  by  integrating 
equation  (2.3): 


Too 

G{A)  =  P{r>A)  =        p{r)dr 


(2.6) 


Figure  3  illustrates  a  Rician  CDF. 
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Cumulative  Distribution  Function  (CDF) 


Usually  the  CDF  is  plotted  in  Rayleigh  paper  which  has  the  especial  character- 
istic that  a  Rayleigh  curve  drawn  in  it  is  a  straight  line.  This  is  useful  for  evaluating  any 
outcome  by  visually  comparing  it  with  the  Rayleigh  line.  Figure  4  illustrates  the  same 
Rician  CDF  seen  in  Figure  3  but  on  probability  paper. 
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Figure  4.      Cumulative  Distribution  Function  (CDF)  on  probability  paper. 


D.     ATMOSPHERIC  EFFECTS  ON  THE  FADING  SIGNAL  COMPONENTS 


The  direct  and  the  difiiise  components  are  both  affected  by  Ionospheric  and 
Tropospheric  effects.  These  effects  include  Faraday  rotation,  group  delay,  absorption, 
dispersion,  refraction  and  scintillation,  all  of  which  result  from  interaction  with  the 
earth's  magnetic  field  and  the  ambient  electron  content  as  the  signal  passes  through  the 
Ionosphere.  Table  1  [Ref  8]  summarizes  some  typical  maximun  values  obtained  for  an 
elevation  angle  of  30  °.  one  way  propagation,  at  two  different  frequencies.  Faraday  ro- 
tation can  be  minimized  using  circular  polarization.  Scintillation  can  be  important  at  low 


latitudes  (9°  N  to  21  °  S)  and  auroral  latitudes  (above  59  °  N),  but  they  are  insignificant 
for  elevation  angles  greather  than  10  °  and  frequencies  below  10  Ghz  [Ref  9]. 


Table   1.     ESTIMATED  MAXIMUM  IONOSPHERIC  EFFECTS  [REF.  8]:     For 

an  elevation  angle  of  30  °  ,  one  way  propagation,  and  a  zenith  electron 
column  of  10'*  electrons'm^ 


Effect 

Frequency 
dependence 

Magnitude 

850  MHz 

1600  MHz 

Faraday  Rotation 

1  '/ 

150° 

42° 

Group  delay 

1    / 

0.35  s 

0.1  s 

Variation  in  direction  of  arrival 

1  /  / 

16  sec  of  arc 

4.7  sec  of  arc 

Refraction 

1  ■/ 

>  50" 

>  14" 

Absorption  (niid-lats) 

\:f 

>  0.014  dB 

>  0.004  dB 

Dispersion 

1   '/ 

0.65  nsec'MHz. 

0.1  nsec  MHz 

Tropospheric  effects  result  from  moisture  in  the  atmosphere  and  tend  to  become 
more  severe  as  frequency  increases.  Table  2  [Ref  10]  summarizes  some  estimated  values 
for  an  elevation  angle  of  30  °,  one  way  propagation,  at  two  different  frequencies. 

The  antenna  pattern  will  influence  the  magnitude  of  the  diffuse  component.  The 
antenna  gain  rolls  olT  below  the  horizon;  therefore,  the  contribution  of  the  diffuse  com- 
ponent is  due  to  scattered  waves  arriving  from  angles  above  the  horizon. 

The  current  treatment  assumes  that  the  diffuse  component  is  8  to  14  dB  (rms)  below 
the  direct  component  [Ref  9]. 
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Table  2.     ESTIMATED  TROPOSPHERIC  ATTENUATION  [REF.  10]:     For  an 

elevation  angle  of  30  °,  and  one  way  propagation 


Effect 

Magnitude  (dB) 

850  MHz 

1600  MHz 

Clear  air  absorption 

3  g  /r?^  (dry) 

0.06 

0.07 

7.5  g'm^  (average) 

0.06 

0.07 

17  g  /?;'  (moist) 

0.06 

0.07 

Cloud  attenuation 

0.5  gw-\  1  Km  thick 

<0.01 

<0.01 

1  g  w-\  2  Km  thick 

<0.01 

<  0.01 

Fog  attenuation 

0.05  g  nv  (average).  0  to  75  Kmt  ht. 

0.05  %'nv  (hea\y).  0  to  150  mt  ht. 

Rain  attenuation 

5  mm  h 

<  0.01 

<  0.01 

25  mm  h 

<  0.1 

<0.1 
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III.     SOFTWARE  FADE  SIMULATOR  FOR  A  MARITIME  MOBILE 

COMMUNICATION  CHANNEL 


A.  INTRODUCTION 

Arredondo  et  al.  [Ref.  11]  describe  a  hardware  simulator  to  simulate  multipath  fad- 
ing for  mobile  radio,  and  Karim  [Ref  12]  describes  a  different  solution  to  implement  this 
simulator  using  software.  This  thesis  follows  both  reference  models,  adding  a  constant 
signal  (direct  component)  to  obtain  Rician  statistics.  It  emphasises  the  doppler  effect 
caused  bv  the  hich  velocitv  of  low  altitude  satellites  and  the  vehicle's  velocity.  Figure  5 
is  a  block  diagram  of  the  suggested  model. 

B.  SIMULATOR  DESCRIPTION 

1.     Rayleigh  Multipath  Fading  Generator. 

The  most  important  part  of  the  simulator  is  the  Rayleigh  multipath  fading  com- 
ponent. Its  implementation  follows  the  model  of  Arredondo  et  al.  [Ref  11]  and  Karim 
[Ref  12].  Two  independent  Gaussian  sources  with  zero  mean  and  equal  variance  are 
generated.  They  are  passed  through  a  shaping  filter  and  finally  added  in  quadrature  to 
obtain  a  signal  with  an  envelope  that  is  Rayleigh  distributed  in  amplitude  and  uniformly 
distributed  in  phase. 

The  generator  is  illustrated  in  Figure  6.  A  key  aspect  in  the  design  of  this  gen- 
erator is  the  shaping  filter. 
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Figure  5. 


Riciaii  fading  simulator  block  diagram 
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A   RAYLEIGH  FADING  SIMULATOR 
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Figure  6.      Rayleigh  fading  generator  block  diagram 

Shaping  Filter:  The  spectral  density  S{a;)  of  the  envelope  of  the  Rayleigh  fading 
signal  received  by  an  omnidireciional  antenna  is  given  by  [Ref  12]  : 


a 


^(<")=-S7n-( 


(JJ  —  (.0 


0   .l~,-0.5 


CO. 


)'] 


\(X)  -  0)^^  <  (Oj 


(3.1) 


where 
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•  (o„  is  the  carrier  frequency. 

•  (x)^  is  the  doppler  shift  due  to  the  velocity  of  the  vehicle  and  the  velocity  of  the 
satellite. 

•  G^  is  the  variance  of  the  two  independent  Gaussian  sequences. 

Refering  to  equation  (3.1),  taking  co  relative  to  oj„  and  noting  that  we  are  only 
interested  in  relative  values  of  the  signal,  the  transfer  function  H(w)  that  has  to  be  ob- 
tained to  shape  the  input  is  given  by: 

H{aj)  =  l\-{^fT'-''  (3.2) 

This  function  is  unbounded  for  co  =  co^,  but  following  the  model  of  Arredondo 
et  al.  [Ref  1 1],  the  actual  transfer  function  has  a  maximum  at  0.95(Dj  and  falls  to  0  with 
a  slope  of -18  dB  ;  octave.  It  is  done  this  way  so  that  the  output  envelope  of  the  simu- 
lator closelv  fits  a  Ravleish  distribution. 

This  transfer  function  will  be  implemented  as  the  transfer  function  H(z)  of  a 
non-recursi\e  digital  filter  with  sampling  frequency  w^  =  2(Oj,  and  the  procedure  is  dis- 
cussed in  Karim  model  [Ref  12]. 

20 

//(2)  =  ^JI>,_,o2-'  (3.3) 

/=0 

where 

•  IV,  is  the  Hammins  window. 
a,  is  the  Fourier  coefficient  of  the  Fourier  series  expansion  of  equation  (3.2). 

z  is  expUcoT)  where  T=  — |-  . 

Figure  7  illustrates  both  the  theoretical  and  actual  transfer  functions. 


• 
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Figure  7.        a)    Spectral    density    of   complex    envelope,    b)    Shaping    filter    freq. 
response.:     [Ref.  11] 

Doppler  Shift  [Ref.  13] 

Doppler  shift  is  caused  by  a  length  change  with  time  in  tlie  propagation  path. 
Let  the  signal  from  the  satellite  be  A^  cos{2nfQt).  This  signal  arrives  to  the  receiver  after 
a  time  delay  g{i),  where  g  is  a  function  describing  the  time  variation  in  path  length. 
Therefore  the  received  signal  is  given  by  Aocos[2nfo{t±g{t))].  Differentiating  the  phase 
we  can  obtain  the  frequency  of  the  received  wave,  that  is: 


1      d  (       Ml) 


(3.4) 


Thus  the  shift  frequency  is: 
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A/=^t^-^  (3.5) 


D{t) 
Since  g{t)  =     ^     where 

•  D{t)  is  the  slant  range,  and 

•  C  is  the  speed  of  light, 


dg{t)       d[D{t)jC\        , 
we  have 


dt  di  C 

where  v  is  the  relative  velocity  of  the  receiver  with  respect  to  the  source.  Due  to  the 

higher  velocity  of  the  satelHte  compared  with  the  velocity  of  the  receiver  on  earth,  the 
latter  becomes  negligible.    For  the  purpose  of  this  thesis,  only  the  radial  component  of 
the  velocity  of  the  satellite  v;  will  be  considered. 
Therefore  equation  (3.5)  becomes: 

A/=+/o-^  (3.5a) 


For  satellites 


d{t)  =  liH  +  Rf  +  R;  -  2R,{H  +  R,)  cos  0]  (3.6) 


where 

•  H  is  the  altitude  of  the  satellite  orbit 

•  R,  is  the  radius  of  earth  =  6378.155  Km 

•  ^  is  a  linear  function  of  time  for  circular  orbits. 

Due  to  the  fact  that  earth  terminals  are  not  located  at  the  center  of  the  earth,, 
there  is  a  maximum  Doppler  shift  at  maximum  slant  range  (satellite  "rise"  and  "set"),  as 
shown  in  Figure  8.  From  the  same  figure  the  following  are  obtained: 

R,  =  R,  +  H  (3.7) 

y  =  sin-'(^^^  (3.8) 

vy  =  V  sin  y  (3.9) 


17 


where 

•  R,  is  the  total  radius  of  the  satellite  orbit. 

•  y  is  the  angle  described  in  Figure  8 

•  \>  is  the  radial  velocity  of  the  satellite,  from  which  the  maximum  value  of  Doppler 
shift  is  obtained. 


Figure  8.      Doppler  shift  from  a  Ion  altitude  satellite 


Therefore,  the  maximum  value  of  Doppler  shift  is: 


^r,max   ^  ^r,max 


fd.max  —  i/ol        Q       )   —    ± 


;. 


(3.10) 


where 
•     y^  is  the  doppler  frequency 
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•  /  is  the  wavelength  of  the  carrier  frequency 

•  y^  is  the  transmitted  carrier  frequency. 

For  a  circular  orbit  satellite,  the  velocity  is  given  by: 


v  =  7X  <^"' 


where 

Kmr' 


• 


M,^  ^(w„,„  +  mj  :^  g{m,,J  ^  3.986013x10^ 


sec^ 


From  equations  (3.9),  (3.10)  and  (3.11)  we  obtain: 


Amax  =  ±fo(-ct)  =  ±fo(—^^—)^±  T^  (3.12) 


"■J  ^  *  ^5 


Rs^ 


2.     Direct  Component  part. 

Let  the  direct  component  of  the  signal  transmitted  by  the  satellite  be 
Aq  cos{o)Qt),  The  signal  at  the  receiver  will  be  affected  by  the  velocity  of  the  satellite,  the 
velocity  of  the  vehicle,  or  both,  causing  a  Doppler  shift.  Davarian  [Ref  14]  suggests  a 
hardware  solution  to  take  care  of  the  Doppler  shift  in  the  direct  component:  this  is 
shown  in  Figure  10.  As  a  result,  the  signal  at  the  receiver  will  be  AqCos[{coq  +  (ji>j)t]  at 
satellite  rise  and  .^o'^^'^K^^o  —  (^j)t]  at  satellite  set.  One  of  these  signals  is  added  to  the 
Rayleish  fadine  sisnal  to  obtain  the  Rician  fadin2  Sisnal. 
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Figure  9.      Block  diagram  for  frequency  sliifting  of  the  direct  patli  [Ref.    14  j:     cOo  is 

the  angular  frequency  of  the  carrier  and  co^  denotes  the  Doppler  shift. 


C.     COMPUTER  PROGRAM 

A  computer  program  was  developed  to  simulate  Rician  fading  statistics  for  a  mari- 
time mobile  satellite  communications  channel.  The  computer  program  is  written  in 
Fortran  and  is  included  as  Appendix  A. 
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D.     ANALYSIS  OF  RESULTS 

The  simulation  of  several  sequences  of  multipath  Rician  fading  was  done  assuming 
a  communication  mobile  maritime  link  with  a  satellite  in  circular  orbit  at  an  altitude  of 
500  nautical  miles.  As  was  stated  earlier  in  this  chapter,  the  speed  of  the  maritime  vehicle 
is  considered  negligible  in  comparison  to  the  speed  of  the  satellite.  As  a  result,  the  vehi- 
cle's speed  is  not  taken  into  account.  Figures  10  through  12  show  the  generated  Rician 
sequences  for  values  of  K  (Rician  factor)  equal  to  0,  5  and  10  dB,  a  frequency  of  869 
MHz  and  a  velocity  of  30  knots.  Figures  13  through  15  illustrate  the  histogram  of  the 
preceeding  sequences.  Figures  16  through  31  show  the  cumulative  distribution  function 
(CDF),  using  probability  paper,  of  several  multipath  sequences  for  the  different  combi- 
nations of  the  following  parameters  K  =  -26,  -10,  0,  10  dB,  frequency  =  869,  1501 
MHz,  and  velocity  =  0,  14357  Kts. 

In  order  to  analize  the  generated  Rician  sequences  we  look  at  tw^o  difTerent  CDF 
graphs.  Comparing  Figures  16  and  28  using  the  same  frequency  and  velocity  for  values 
of  K  =  -26  dB  and  K  =  10  dB,  it  can  be  noticed  that  for  the  former  value  of  K  the  re- 
ceived signal  is  greater  than  the  abscissa  for  a  longer  period  of  time.  This  means  that  for 
a  larger  Rice  factor  (K)  we  have  more  fading  and  therefore  less  strength  in  the  received 
signal  at  the  antenna. 

Additionally,  the  difference  in  the  CDF  can  be  noticed  between  the  curves  with  0 
velocity  (0  doppler)  and  the  curves  with  a  velocity  of  14357  Kts  (max.  doppler);  as  was 
expected,  the  difference  is  larger  at  higher  frequencies. 
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Figure  10.      Rician  fading  sequence  for  K  =  0  dB,  f  =  869  MHz,  V  =  30  Kts. 
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Figure   11.      Rician  fading  sequence  for  K  =  5  dB,  f  =  869  MHz,  V  =  30  Kts. 
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K  =   10  dB,  f  =  869  MHz,  V,  =  30  Kts 
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Figure  12.      Rician  fading  sequence  for  K  =   10  dB,  f  =  869  MHz,  V  =  30  Kts. 
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Figure   13. 


Histogram  of  a  Rician  fading  sequence  for  K  =  0  dB,  f  =  869  MHz,  V 

=  30  Kts. 
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Figure   14.      Histogram  of  a  Rician  fading  sequence  for  K  =  5  dB,  f  =  869  MHz,  V 

=  30  Kts. 
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K  =   10  dB,  f  =  869  MHz,  V,  =  30  Kts 
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Figure  15.      Histogram  of  a  Rician  fading  sequence  for  K  =    10  dB,  f 

V  =  30  Kts. 
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Figure   16.      CDF  of  a  Rician  fading  sequence  for  K  =  -26  dB,  f  =  869  MHz,  V 
OKts. 
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Figure  17.      CDF  of  a  Rician  fading  sequence  for  K  =  -26  dB,  f  =  869  MHz,  V 
14357  Kts. 
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Figure   18.      CDF  of  a  Rician  fading  sequence  for  K  =  -26  dB,  f  =    1501  MHz,  V 
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Figure   19.      CDF  of  a  Rician  fading  sequence  for  K  =  -26  dB,  f  =    1501  MHz,  V 
=  14357  Kts. 
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Figure  20.      CDF  of  a  Rician  fading  sequence  for  K  =  -10  dB,  f  =  869  MHz,  V 
OKts. 
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Figure  21.      CDF  of  a  Rician  fading  sequence  for  K  =  -10  dB,  f  =  869  MHz,  V 
14357  Kts. 
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Figure  22.      CDF  of  a  Rician  fading  sequence  for  K  =  -10  dB,  f  =    1501  MHz,  V 
=  0  Kts. 
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Figure  23.      CDF  of  a  Rician  fading  sequence  for  K  =  -10  dB,  f  =    1501  MHz,  V 
=  14357  Kts. 
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Figure  24.      CDF  of  a  Rician  fading  sequence  for  K  =  0  dB,  f  =  869  MHz,  V  =  0 
Kts. 


36 


99.9 


<      99 

V) 
CO 

< 

A      75 


50 


^      25 


< 
z 
o 


a 
> 

UI 

u 


10 
5 

1 


3     0.1 

I- 

H-    0.01 

z 

yo.001 


K  =  0  dB.  f  =  869  MHz,  V,  =   U357  Kts 

CDF 


\                      •                      ■        /  '          :                      ' 

\     \    y    \     I 

:                      j 

''■■'■    J  •■     '■■ 

I                1               j 

.i^f... i 1 i 1. 

I \,^r^ 

y^f          \            T            i            ] 

1,.^1 

J^-          \ 



\                        \                        \                        '.                        \ 

-15  -10  -5  0  5 

RECEIVED  SIGNAL  RELATIVE  TO  LINE  OF  SIGHT  (dB) 


Figure  25.      CDF  of  a  Rician  fading  sequence  for  K  =  0  dB,  f  =  869  MHz,  V 
14357  Kts. 
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Figure  26.      CDF  of  a  Rician  fading  sequence  for  K  =  0  dB,  f  =   1501  MHz,  V  = 
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Figure  27.      CDF  of  a  Rician  fading  sequence  for  K  =  0  dB,  f  =    1501  MHz,  V 
14357  Kts. 
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Figure  28.      CDF  of  a  Rician  fading  sequence  for  K  =   10  dB,  f  =  869  MHz,  V 
OKts. 
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Figure  29.      CDF  of  a  Rician  fading  sequence  for  K  =   10  dB,  f  =  869  MHz,  V 
14357  Kts. 
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Figure  30.      CDF  of  a  Rician  fading  sequence  for  K  =   10  dB,  f  =   1501  MHz,  V 
OKts. 
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Figure  31.      CDF  of  a  Rician  fading  sequence  for  K  =   10  dB,  f  =   1501  MHz,  V 
14357  Kts. 
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IV.     CAPACITY  OF  UNSLOTTED  ALOHA  MARITIME  MOBILE 

SATELLITE  SYSTEMS 

A.  INTRODUCTION 

Aloha  is  a  multiple  access  method  commonly  used  in  data  communications.  This 
chapter  analyzes  the  throughput  of  unslotted  Aloha  with  power  capture  in  a  Rician 
fading  environment. 

B.  UNSLOTTED  ALOHA  THROUGHPUT 

In  the  Aloha  system,  a  satellite  channel  with  a  capacity  of  R  bits  per  second  is 
shared  by  a  large  population  of  users.  Each  user's  station  randomly  transmits  packets 
at  the  channel  bit  rate  R  whenever  its  buffer  contains  one  packet.  When  packets  from 
different  stations  overlap  (packet  collision),  the  transmission  error  can  be  detected  at  the 
transmitting  stations  on  the  downlink.  The  stations  then  retransmit  the  packets  until 
they  are  free  from  overlap.  To  avoid  repeated  collisions,  the  interval  of  packet  trans- 
mission is  randomized  for  each  station  [Ref  15]. 

In  unslotted  Aloha  there  are  two  types  of  capture  effect:  created  and  natural.  Cre- 
ated capture  effects  occur  when  random  power  levels  are  used.  In  this  case,  the  packet 
with  the  higher  power  level  may  capture  the  receiver  when  its  signal-to-interference  ratio 
exceeds  a  preset  threshold.  Natural  capture  eflects  occur  in  mobile  radio  communi- 
cations when  varying  distances  between  mobile  nodes  and  base  stations,  as  well  as 
channel  fading,  cause  var\'ing  signal  powers  at  the  base  station's  receiver  [Ref    16  ]. 

As  stated  in  the  paper  by  Borchardt  et  al.  [Ref  17]  about  unslotted  Aloha,  let  us 
assume  an  infinite  number  of  uncoordinated  users  communicating  to  a  base  station  us- 
ing fixed  length  packets.  The  channel  traffic  is  approximated  by  a  Poisson  proccess  with 
parameter  G  packets  per  packet  length.  Thus  the  probability  that  a  tagged  packet  is 
overlapped  by  n  other  packets  during  its  transmission  is  the  arrival  probability  of  n 
packets  during  the  interval  {i^—  T,Iq+  T)  ,  where  i^  is  the  transmission  instant  of  the 
tagged  packet  and  T  is  the  packet  transmission  time.  This  arrival  probability  is 

Pr{«}=-^^exp(-2G)  (4.1) 

n. 
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According  to  Borchardt  et  al.  [Ref.  17],  consider  a  receiver  that  captures  a  tagged 
packet  in  the  presence  of  n  interfering  packets,  if  the  tagged  packet  power  sufficiently 
exceeds  the  joint  interference  power.  These  interferers  form  a  stochastic  proccess  which 
consists  of  early  interferers,  whose  transmission  began  in  the  interval  (/g  —  1,  ^o).  arid  late 
interferers  which  enter  the  system  in  the  interval  (ro/o  +  1).  With  n  interferers  there  are 
2"  realizations  of  early  and  late  interferers.  C{n)  is  the  number  of  realizations  in  which 
the  maximum  number  of  interferers  does  not  exceed  N  (the  maximum  number  of 
interferers  that  can  be  tolerated  at  any  instant  of  time  during  the  tagged  packet  trans- 
mission) and  is  given  by  [Ref  18] 

C[n)=    ^  Cj{n)       ,  [jc]  =  integer>;c  (4.2) 

y=[f] 


where 


n\  (2/-A7+  1)^ 
7+1 


,J>nl2 
,  otherwise 


where  C,{n)  is  the  number  of  realizations  in  which  the  maximum  number  of  interferers 
equals/ 

The  capture  effect  in  a  random  access  channel  occurs  naturally  when  fading  reduces 
the  joint  power  of  interfering  packets,  as  reported  in  [Ref  19]  for  a  slotted  Aloha  channel 
in  a  uniform  Rayleigh  fading  environment.  This  thesis  extends  these  results  to  unslotted 
Aloha  in  a  Rician  fading  environment. 

It  is  assumed  that  a  tagged  packet  will  capture  the  receiver  if  its  signal-to- 
interference  ratio  exceeds  a  threshold  y^. 

The  probability  density  function  of  a  packet's  power  in  a  Rician  fading  channel  is 

fpip)  =  —  k\^ a )  ^^P(^ a )  (4-3) 

The  same  expression  for  a  normalized  mean  power  a  =  1  is 

fpip)  =  exp(  -^oVo(2^ov>  )  exp(  -p)  (4.4) 
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Thus  given  n  interfering  packets  during  the  transmission  of  a  tagged  packet,  the 
density  function  of  the  power  of  n  interfering  packets  is  [Ref  17] 


/y(yl«)  =  -3r   £  Cj{n){fp{p)y®  (4.5) 

Substituting  equation  (4.4)  into  equation  (4.5)  we  get 


n 

fy<y\")  =  ^r   Z  ^/"^^  ^^P(  -A'M{2A,jy)  exp(  -y)y®  (4.6) 

where  J®  denotes  the  j-fold  convolution. 
Appendix  B  illustrates  the  derivation  of 

[  exp(  -  AS)Io{2AQ^y)  exp(  -y)J®    for  J  =  1,2,3,4,5  and  6 

Y 

The  densitv  function  of  the  power  of  the  sisnal-to-interference  ratio  Z  =  -^  is  ob- 

tained  bv 


/z(^H«)  = 


yfx<y^Vy<y\n)dy  (4.7) 


Thus 


/z(-l«)  =  — 


n 


y  exp(  -  AQ)lQ{2AoJyz  )  exp(  -yzV  2_j  Q(^)l^  ^^P  "  ^o^(2^ov'7)- 

l/=[^/2]  (4.8) 


•  exp(  -y)y®} 

The  capture  probabiUty  is  given  by 

Pr{capture\n,  n>\]  =  1  -  /^,{0<Z<yol«,  n>l]  (4.9) 

Therefore  we  have 

fvo 
Pr{capture\n,  n>l]  =  1  -      fzi^W^  (4T0) 
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Appendix  B  illustrates  the  derivation  of  the  probabilities  of  capture  for  n  =  1,2,3,4,5 
and  6. 

Channel  throughput  S  is  defined  by 


S  =  GPr{capture]  =  G/  Pr{capture\n]Pr{n] 


(4.11) 


n=Q 


Replacing  equations  (4.1)  and  (4.9)  in  equation  (4.11),  the  throughput  of  unslotted 
Aloha  in  a  Rician  fading  channel  is  given  by 


-"E^ 


S=G 


exp(  —2G)Pr{capture\n] 


(4.12) 


/7=0 


and  given  that  Pr[capture\0]  =  1  it  reduces  to 


G  exp(  -2G) 


oo 

Z2G" 
— ^  Pr[capture\n] 


H=l 


(4.13) 


Figures  32.  33  and  34  illustrate  the  throughput  for  several  capture  thresholds  y^  us- 
ing y4o  =  0.1,  Aq  =  0.5  and  A^  =  1  respectively.  Finally,  Figure  35  illustrates  the  through- 
put for  the  same  capture  thresholds  >o.  but  letting  parameter  Aq  =  Q  (Rayleigh  fading). 
As  expected  the  same  throughput  of  unslotted  Aloha  in  a  Rayleigh  fading  environment 
[Ref  17]  were  obtained. 
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V.     CONCLUSIONS 

This  thesis  can  be  divided  into  two  main  parts.  The  first  part  describes  the  imple- 
mentation of  a  software  simulator  for  multipath  fading  in  a  maritime  mobile  satellite 
channel.  The  second  part  analyzes  the  throughput  for  an  unslotted  Aloha  channel  in  a 
Rician  fading  environment. 

The  fading  simulator  generates  Rician  fading  sequences  for  any  value  of  frequency, 
velocity  and  Rice  Factor  (K).  Fades  can  be  simulated  for  communication  links  using  a 
geostationar\'  satellite  or  low  altitude  satellites. 

Unfortunately,  it  was  not  possible  to  obtain  data  from  experiments  done  in  several 
institutions  like  JPL  and  CRC  in  order  to  compare  results  and  make  necessarv'  adjust- 
ments to  the  simulator. 

For  the  second  part  of  the  thesis,  a  set  of  equations  are  derived  to  compute  the 
throughput  of  unslotted  Aloha  in  a  Rician  fading  channel.  The  obtained  equations  are 
tested  for  several  values  of  the  parameters  A^  and  y^.  Finally,  a  comparison  is  made  with 
results  obtained  for  a  Rayleigh  fading  channel  [Ref  17]  by  setting  the  parameter  A^^  equal 
to  zero,  and  similar  results  are  obtained.  It  is  noted  that  the  throughput  of  a  Rician 
fading  Aloha  channel  depends  on  the  direct  component  power  A^.  As  Al  increases,  the 
joint  interference  power  also  increases  and  consequently  the  throughput  is  reduced. 
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APPENDIX  A.      SOFTWARE  RICIAN  SIMULATOR 

This  program  generates  a  Rician  fading  sequence  given  variance,  frequency  and  ve- 
locitv. 


C      yf*,v*-c  VARIABLE  DEFINITIONS  ***** 

C      A   =   LOVER  LIMIT 

C      B   =   UPPER  LIMIT 

C      FCT  =   FUNCTION 

C      S   =   PARTIAL  SUM  VALUE 

C      H   =   INTERVAL  LENGTH 

C      AC  =   SIMPSON'S  RULE  VALUE 

C      N   =   NUMBER  OF  EVEN  SUB INTERVAL 

C     NN2,NN,I,CC    =  COUNTERS 

C      X,H,D         =  COMPONENTS  OF  THE  RULES 

C     VARIABLE  DECLARATIONS 

COMPLEX  RF(0: 1024), R(0:  1024),DSV(0:  1024) 

REAL  SI(0:  1024), SQ(0:  1024) , V,ME ,G( 0:  500) ,F( 0: 500) , Y,W,CONST( 0:  1000 
&),WO,DS 

REAL'^^S  FR,WS,WL,A,B,FCT,S,X,H,D,AC(21)  ,PI,Z,RGF(0:  1000)  ,AX(0:  1000) 
&,IGF(0:  1000),MAG(0:  1024) ,TS 

INTEGER  I,NN2,N,NN,NS 
C     MAIN  PROGRAM 
C     Nl: NUMBER  OF  SAMPLES  OF  THE  GAUSSIAN  RANDOM  GENERATOR 

N 1=5 00 
C     N2: NUMBER  OF  COEFFICIENTS  CALCULATED  FOR  THE  SHAPING  FILTER 

N2=21 

NP  =  Nl  +  N2 
C 
C     GAUSSIAN  RANDOM  GENERATOR  1  WITH  VARIANCE  1,  MEAN  0 

ISEED=23 

VRITE(--'S^'0  'ENTER  VARIANCE' 

RE  AD  (•'-,-'■)  V 

ME=0.  0 

DO  10  1=0, Nl 

CALL  GAUSS  ( ISEED,V,ME,W) 
G(I)=W 

10  CONTINUE 

C     GAUSSIAN  RANDOM  GENERATOR  2  WITH  VARIANCE  1,  MEAN  0 
DO  11  L=0,N1 

CALL  GAUSS  ( ISEED, V,ME , Y) 
F(L)=Y 

11  CONTINUE 

C     WRITE(8,*)  (G(I),I=0,N1) 

C     WRITE(4,*)  (F(I),I=0,N1) 

C 

C     CALCULATION  OF  SHAPING  FILTER  COEFFICIENTS 

C 

c 

C     TO  CALCULATE  A  DEFINITE  INTEGRAL  FOR  A  FUNCTION  USING  SIMPSON'S 
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C     RULE. 
C 

NS=21 

PI=DATAN(1D00)''^4 
C=2. 997925E+08 

WRITE( '•-/-)  'ENTER  OPERATING  FREQ' 
READ(*,*)FR 
WO  =  2  *  PI  ''■-  FR 

WRITE(*,*)  'ENTER  SPEED  IN  KNOTS' 
READC^.'-OWM 
IF  (WM  .EQ.  0)  THEN 
DO  12  K=0,N1 

DS  =  COS(WO-VK) 

RF(K)  =  CMPLX(G(K),F(K)) 

DSV(K)  =  CMPLX(DS,0. 0) 

R(K)  =  DSV(K)  +  RF(K) 

MAG(K)  =  CABS(R(K)) 

WRITE (4,''-)  MAG(K) 

12  CONTINUE 
ELSE 
WL=C/FR 

WRITEC'S'-O  'WAVELENGTH  =  '  ,WL 
WM=2-'-PI''^(WM'n852/3600)/WL 
WS=WM/2/PI 

WRITE('>,^'0  'DOPPLER  FREQ  =  '  ,WS 
B=0.9  7''-WM 
WS=2-''^WM 
TS=2^'^PI/WS 

WRITE (--'-,'■-)  'ENTER  (EVEN)  N  FOR  SIMPSON  INTEGRATION' 
RE  AD  (--••-,  ^'ON 
C 

A=0 

H=(B-A)/N 

D=H/3D00 

WRITE( ''-,*)  '^"''"'"'^''^COMPUTING:  PLEASE  WAIT*****'^*' 

DO  13  K=0,  (NS-l)/2 

S=FCT(WL,K,WS) 

DO  14  1=1, N-1,  2 
WL=A+I''^H 
S=S+4.  O^^FCT(WL,K,WS) 

14  CONTINUE 

DO  15  1=2,  N-2,  2 
WL=A+I--''H 
S=S+2.0'^FCT(WL,K,WS) 

15  CONTINUE 
S=S+FCT(B,K,WS) 
AC(K+1+(NS-1)/2)=S''^D 

13  CONTINUE 

DO  20  I=l,(NS-l)/2 
20    AC(I)=AC(NS+1-I) 

DO  21  J=0,(NS-1) 
AX(J)=AC(J+1) 
21     CONTINUE 

DO  25  N=0,NS-1 
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AX(N)=AX(N)*(0.54-0.46*DC0S(2*PI*N/(NS-1))) 
25    CONTINUE 
C 

C     CONVOLUTION  OF  GAUSSIAN  RANDOM  VECTOR  1  AND  SHAPING  FILTER  1 
DO  26  N=N2,N1+N2 
AX(N)=0.  OD+00 
26   CONTINUE 
C 

CALL  CONVOL  (N1,N2,G,AX,RGF) 
C 
C     CONVOLUTION  OF  GAUSSIAN  RANDOM  VECTOR  2  AND  SHAPING  FILTER  2 

CALL  CONVOL  (N1,N2,F,AX, IGF) 
C 

C     APLYING  DOPPLER  EFFECT  TO  DIRECT  SIGNAL 
DO  30  K  =  0,NP 
DS  =  COS  ((WO  +  WS)''^K'^TS) 
C 

RF(K)  =  CMPLX(RGF(K),IGF(K)) 
C     ADDING  RAYLEIGH  SIGNAL  TO  A  CONSTANT  TO  OBTAIN  RICEAN 
DSV(K)  =  CMPLX(DS,0.  0) 
R(K)  =  DSV(K)  +  RF(K) 
WRITE(8,*)  R(K) 
MAG(K)  =  CABS(R(K)) 
WRITE(4,''-)  MAG(K) 
30   CONTINUE 
ENDIF 
STOP 
END 
C 

REAL-'^8  FUNCTION  FCT(WL,K,WS) 
REAL-8  WL,WS,PI,D 
INTEGER  K 
PI=DATAN(1D00)-'^4D00 
D=2D00--'-PI-'>K''-WL/WS 

FCT=2D00/WS-"-(lD00-(2D00-'^WL/WS)''"V2)*-V(.0,  25)''^DC0S(D) 
RETURN 
END 
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APPENDIX  B.     DERIVATION  OF  FORMULAS  FOR  UNSLOTTED 

ALOHA  THROUGHPUT 


A.     J-FOLD  CONVOLUTION 

Given  the  expression 


[e-^°/o(2^o.(7)^""] 


{BA) 


in  order  to  compute  the  j-fold  convolution  for  j  =   1,2,  3,  4,  5  and  6,  we  evaluate  the 
characteristic  function  and  obtain 


Oy{u) 


.-^%{2^ov5^)^~V"^^ 


(5.2) 


Replacing  the  modified  Bessel  function  of  order  0  by  its  series  approxir^.ation  to  the 
fourth  power  we  have 


CD,{^)  = 


-a; 


e-yei-y  +  e-'Alye^''  +  e''  -^  j'V"^  +  e''  ^>-V"^ 


—y      '^0  4  JuX 

:>ib  -^ 


{B.3) 


dv 


This  is  equivalent  to 


Oy(«)  =  e 


-a; 


■CX) 

e 

6      /^oo 


-0-M. 


'■"^%  +  ^o'j 


,4     Too 
°°       -{\-Ju)y,     ,     ^0 

ye  "'dy  +  — r- 


2  -(l-yw)^  - 
^  ^  ^    •^     dy 


+ 


'0 


36 


j;  e  '    ^  '■dy  +  -T= 


8    r*oo 
576 


/,-(i-».^^ 


(5.4) 


After  evaluating  the  integrals  we  obtain 
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cDy<u)  =  e 


-4 


+ 


A 


1 

(1-7") 


24 


+  ^ 


2 


1 


+ 


»0 


+ 


'  {\-Juf    '     4     (l-j^)^    ■    36    ^x-juf 


576  (i-j^r 


{B.5) 


Now,  the  j-fold  convolution  in  time  domain  is  equivalent  to  the  jth  power  in  the 

frequency  domain,  that  is 


m-)r'^i<^M)1 


Evaluating  the  different  jth  powers  and  taking  the  inverse  transform  to  each  expression 
we  obtain  the  desired  j-fold  convolutions. 
For  j  =  1 


®  _     -A't 


(fyW  =  e 


j4  ,6 

e-'  +  lA^ve-'  +  0.5  ^y'e-^  +  0.1666-^/^"^ 


a...-y 


{B.6) 


Forj  =  2 


Ar 


A 


-y    1    -»    "'^'       2    ->>    ,    /N    '^0      3    -y    ,    1   -3  -.  ->  -5       0      4    -y 


ye  '  + 

^3 


2,  -/.-^  +  2  -rr/^"-'  +  1-3333  -^ye 


iO.,.5.->' 


+  0.5833  -rry  e~'  +  0. 1666  ^^  ve~'  +  0.0277  — ^  v^ 

•»l«^  /.!»'  /!»' 


,10 

0        6    -y 


...   ^0 


12 


,.'^->' 


o: 


6! 


7! 


(5.7) 


For  j  =  3 
(Ay(v))'® 


-3-4r 


1        2    -><    ,    -,    ^0       3    -V    ,     .   c    ^0       4    - V    ,     .   r    "^^ 

— >•  e  "^  +  3— ^  e  '  +4-5— >•  ^  ^  +  4.5  — 


.8  .IOn 

+  3.25— >-g  ^+1.7^  — 


>-'.-^  + 0.7083-^  A"' 


v4 


12 


8! 


(5.8) 


+  0.2083^>V^  +  0.0416^yV^  +  0.0046^>.-''.-^ 
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For  j  =  4 


4! 


4 
5! 


j^/e-'  +  4  ^  A"^  +  8  -^>-^.-^  +  10.6666  ^/."^ 


:1c 

6! 


+  9.6666  ~:^^  e  ^  +  S  y  e  ^  +  4  —^y  ^     +  2.3333 


^0         10    -V 
A6  .18  .20  (^-9) 

+  0.8958  -^j'".-^  +  0.2685  ^>-' V^  +  0.6018  ^j;^^."^ 

+  0.0092  4V  v''^"'  +  0.00077  4V>-''^"' 
14!  -  1:)!  -^ 


For  j  =  5 


.^^5®  -  ^-S-^c 


(/'yl^'))-^  =  ^ 


1         4    - 


^0^ 


"^    - 


0_    .6    -y 


A 


0_     1    -y 


-fye-'  +  5  -f-ye'  +  12.5  "^  A""  +  20.8333  -yy'e 


A^  ^'°1  A^^ 

+  2^.8jj3  •——  V  e  ■"  +  25.1666---^ — >•  e  -^  +  19.8611    ,^,    v'    ^ 
8!  "  9!  10! 

+  12.91664V^"^"'  + 6-9791  ^j^'V^  +  2.6157^>-'V^     (^-lO) 


A^^  A^^  A^'^ 

+  1.1655--p>;    ^  -"  + 0.353-777- J'    ^  ^  +  0.0848— -j- J'    e^ 

^26  ^28  ^30 

+  0.0254  -^y^\-^  +  0.00096  -r^rJ''^^"^  +  0.000128  -r^y^^^e'^ 
17!  ^  18!  *^  19!  *^ 
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For  j  =  6 


y. 


6! 


1        5    -V    ,    ^    -'^0       6    -y    ,    ,  o    ^0       7    -y    ,    ,^    ^0       8    -v 


7! 


lOx 


Ar\        o    _„  _     _    /iQ    j        in    _„  /j 


8! 

12 


+  51>.15—ye^  +  (>Z.5—^y    £  "^  +  61.4166 -^>-    e  "^ 

+  49.6666  ^y'^e-'  +  34.0208  ^J^''^"^  +  19.8842  ^y''e-' 


,20 


1 22 


,24 


4  A  A' 

,    n  r>A/fn        0  1 5    - y    ,     -,    "^0  16    -y    ,    r.  n  i  o /i        0  H    -> 

+  9.0069 -^j;    ^  ~16r-^"    ^     +^-91^^~pr>'    ^ 

.26  .28  .30 

+  0.3322 -Yg7->-    e  ■'  +  0.1013-|^j;    e^  +  0.0181^^y    e  "^ 

^32  ^34 

+  0.0036  -z^y^^e'^  +  0.000385  ^V/^^'^ 


(5.11) 


B.     PROBABILITY  OF  CAPTURE 

First,  evaluating  equation  (4.8). /^(z  |  «)  for  n  =  1  2,  3,  4,  5  and  6,  the  following  ex- 
pressions are  obtained 

For  n  =   1 


/z(z|l)=Q(l).-^'^'^ 

2 


(1  +  ^)' 


g+1 


+  A'z 


y.K-^ 


^  ■' '  (1+^^) 


gig-  1)        ,      ,4 

rrr  +  ^o- 


g(^-l)(^-2) 


-,(Z?.12) 


^>:^-,^^-^^^4tt^+447->a- 


4    /  /  g-2 

^=3 


(l+2f 


g+1 


3S4^^-^^<^- 


(^-4)!(l+2) 


,5+1 
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For  n  =  2 


/2(7|2)  =  e-'^°^|c,(2) 


4 


.4f'^^^T.-) 


'  +  Alz 


^(^-1) 


5=2 


(1+2) 


,g+i  ■^'^o 


4Z-. 


g(^-l)(g-2)    ,      .  ^3 


g=3 


""-'   (1+zr^ 


+  ^0   36 


Zv 


^: 


+  Q(2)^ 


-/^r 


g=4 
9 


'  te-4)!(l+zr' 


(5.13) 


+  ^0 


4    Z 


4    • 


10 


h{h-  l)(/7-2) 


11 


Z^^''-^  ~Tr~T;^Tr^ + ^^  "k  Z^'^ 


,1 


h=A 


(1+2) 


,.3   '-^/z-4)!(l+z)^^> 


For  n  =  3 


/z(z13)  =  .-^^°^|q(3) 


8 


I 


.=.  *'(1+^)" 


h(h-l) 


1     -.  -/,+  !    +^0' 


h=2 


10 


"tZ^''-^ 


M/j-l)(/j-2) 


II 


/,=4  (1+2) 


+  Ao  36 


■i6"Z.^''-3 


(1+2)' 


,t^   ■•"'  (/2-4)!(l+z)''^^ 


+  C3(3)^ 


-A^ 


■   12 


13 


^  +  ^n'z^^;_, 


/(/-I)       ^,4z' 


(5.14) 


.3    "(1+2)-    ^^"^^^'-'(1+2)-^^°     4 


14 


■Z^i-2 


/(/-!)(/- 2)    ,      .  ,3 


/=5  (1+2) 


15 


*^°  36^-"  (,_,),(,  ^^y., 
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For  n  =  4 


/z(z14)  =  ^-^'^°-^|q(4) 


t.^hh—^  +  A', 


h=2 


^''(1+2^ 


h{h-\)     ,      . 

.=3  -^d.zr''^- 


10 


11 


•  X  Z^^''-2       „   .  .,/>+i       +  ^0  "36"  Z^^''- 


//! 


ft=4  (1+2) 


.1^    ■■  '  (/^-4)!(l+z)''+' 


14 


-/4n 


12 


13 


/(/-I)        _ 

/+1    '  '-o^/_,'V-i  ;,   ,  _^ft+i  ^'^o   4  ■ 


u  {\+zr     t^    (1+2) 


+  A 


4    Z 


15 


^K,_, 


1=5 


(1+2)' 


Z^/-3 


/! 


/(/-!)(/ -2)         ,  ,3 


/+i 


(5.15) 


+  Q{A)e 


-2Af 


17 


§''^777;^^  ''"''S"'^-  oT^  " ''°  "4" 


■'■^■-"     +A' 


18 


'l^"^^-^      (1+2/-      ^^°36M-3^._ 


(/■-4)i(i  +  2y 


+  1 


For  n  =  5 


Mz\5)  =  e-'''^-^lc,i5} 


12 


13 


U       (l+z/ 


^^E^/-i 


/(/-I) 


/'    ,.  /+i    -^^O-^'V-l    ^^    ,      ,/+i    +^0- 


/=4 


14 


15 


--„     '<'-■)('-.,, j4^,„ 


'h^-'  (1 


(1+2)' 

/! 


+  z 


J+i 


t^    •  '  (/-4)!(l+z/^' 


+  Qi5)e' 


16 


17 


18 


19 


0^/  /v-i  ..      :/+i  +  ^0  4  • 


V.     .•(/•- l)(/--2)e^y        ^_^ 

^S^  ''(/■-4)!(i  +  2y- 


M--2-    (!+,/+> 


7=6 


+  C5(5). 


-2/ln 


20 


E 


A^^^w 


1 


21 


(5.16) 


m(m  —  1) 


„=5  ■"■■■  (1 + ^r '  ^ ''°\4^e''"-' (I + ^-r '  ^  ^° '-^ 


22 

.y 


23 


K. 


m{m  —  l)(m  —  2)  ^  z^ 


m=i 


m-2  (i^^)-+l 


Z^^'"-3 


m! 


ffi=8 


(m-4)!(l  +zy 
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For  n  =  6 


/z(2|6)  =  ^-^-'«-^|c3(6) 


p'l^'''^'" 


1 J 


/(/-I) 


+  A 


/=4  (I+2) 


14 


/(/-!)(/ -2)  ,  ,3 

4   /l^"'-2      n^,V+i      "^     °  36 


15 


Z^'-= 


/=5 


+  C4(6)^~^o 


(1  +  2)' 

16 


Z'^'-3 


/> 


t?   ■■'  {l-AW+zf 


y=4  I  i   -^  2/  y=5 

18  .,.         .,,.        ^,  ,      19 


7(/--l) 


+  /i 


4  z 


/ 


+  C^[6)e 


-1A-, 


20 


7=7      (/■-4)'(i+2y 

21 


+1 


(R17) 


E^' 


m=:> 


"(1+zr' 


+  ^0^ 


Z^- 


m{m  —  1) 


,4  z 


m=6 


1    ,.     ,      .^+1    +^04- 


(l+z)' 


TT 


m(m  —  l)(m  —  2) 


23 


Z  rriyrri  —    i  jyrri  —  ^,       ^       ^g 


ml 


m+\ 


m=l 


+  Qi6)e-'^^ 


(1+z) 

23 


^^0     36     Z-J^""-^    /...         ,M/i     ,     _^"l+l 


1 


.,=8  (m-4)!(l+z) 

24 


+  AnZ 


i.=6  (!+-")  u=7 


3^7    A'j,_i 


^(^-   1)  4_Zl 


u[u  -  1){U  —  2)  6    z      V*   ^ 


u= 


?/: 


where  the  C,(«)  values  are  defined  in  equation  (4.3),  and  K^,  K^,  K„  Kj,  K„,  SiudK^  are  given 
in  TabJe  3. 
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Table  3.     CONSTANT   VALUES    USED    IN   EQUATIONS   (B.12)   THROUGH 
(B.23) 


a; 

a; 

A, 

a; 

A,„ 

A„ 

1 

1 

1 

1 

1 

1 

1 

2 

2Al 

2Al 

^Al 

^Al 

5^0^ 

(^Al 

3 

O.SAt 

1A\ 

4.5^S 

"^Al 

12.5^g 

18/ig 

4 

0.1666^* 

1.3333.4o' 

^.SAl 

10.1666^^ 

20.8333^* 

36/i« 

5 

0.5833.-^1 

■h.lbAl 

9.6666/}^ 

25.8333/^0* 

53.75.^0' 

6 

0.1666^^0 

l.lSAf 

%Af 

25.1666/^^° 

63.5/4^0 

7 

0.0277^^^ 

0.7083.^'^ 

4/^^^ 

19.8611/^^2 

61.4166/^'^ 

8 

0.2083.4,r 

2.3333/i^^ 

12.9166.^^* 

49.6666.^r 

9 

0.041 6.^  >^ 

0.8958.^^^ 

6.9791.4^^ 

34.0208/^^,^ 

10 

0.0046.Jri' 

0.2685/i^^ 

2.6157.-ii» 

19.8842///,' 

11 

0.0602.-^5'" 

1.1655.^5° 

9.0069.-l5'° 

12 

0.0092.4o^-^ 

0.353.^5'-^ 

O.lAf 

13 

0.00077/lo^^ 

0.0848.^o'' 

0.9284.^,3^ 

14 

0.0254.^5'^ 

0.3322.45" 

.15 

0.00096.-lf 

0.1013.-/,^ 

16 

0.0001 3.-^0^° 

0.018^5^° 

17 

0.0036.-/,^^ 

IS 

0.00038.-/,^^ 
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Now,  the  probabilities  of  capture  are  evaluated  by  replacing  equations  (B.12) 
through  (B.17)  into  equation  {4.10),  and  we  get 


For  n  =  1 


Pr{capiure\l}  =  1  -  ]  -^  e'^'^'c^il) 


Zv(^i7n^-'j)-^ogv. 


■^(^-1) 


-1 


[ 


^^-1)^(1  + -/or' 


(---)      (l+7o)' 


4      6 


'^  °=3 


-1  r  1  '> 


is-i)\i  +  y,r^    '^ 


(5.  IS) 


3 


+ 


-^j!  V  (3-^)    (1  +  yo)^       (3-^) 


^S-^)      (l+Vof 


0—2 


-1]- 


(3-^)     (^-1)        (l+yj 


o— 1 


1]  + 


1 


(3-.-)    S 


•[ 


1 


(i  +  7or 


r-i] 
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For  n  =  2 


Pr{capture\2}  =  1  — 


■^(^-1) 


rr^"'^°Q(2) 


Zv(^i7rT^-n)..o4v 


(^-^)     (i  +  yo) 


^  ;ir-i]  +  i[ 


^    (i  +  yo) 

2 


'    r-ilV^Zv 


•^(^-i)te-<(i^i777^-'i+17^i 


1 


-^[ 


^^^(1  +  yof 


— i — [ 


(^-1)      (1+yof 
1 


4    /  ."g-2- 

^=3 


-1] 


36  /-.   s-3  ^g_4y_   \{3-g)    [l  +  y,f  ^  (3-g)  • 


(^-2)      (1+yof 


-2 


1] 


■[ 


I 


(1  +  yor 


-11 


+ 


1      -3a'; 


(3-.e)   (.^-1)^(1  +  ,^) 

8 


L,     V.-,  -i]  +  7r^i 


(3-.-)   ^ 


C2(2) 


z 

/j=2 


^a/J 


-1 


1 


^'     (1  +  yo/ 


1]    +^c 


(5.19) 


•Z^-^'^^-^^TA^f 


/3=3 

10 


(''-')    (1  +  yo) 


/i-l 


"Hi77t:^-")-4 


(l  +  zo) 


.^AV2M/^-l)(/^-2)(-— ^-[ 


/7=4 


1 


3 


6     11 


-1]  + 


[ 


1 


[^^-I]).-Ag^^,_3_^(_ 


(^-1)     (1  +  yo) 
1 


h-\ 


-1] 


3 

yo 


''^  (i  +  yo)' 


+  TT 


1 


[ 


1 


4- 


(3-/0    (A-2)^(i+,^)^-2 
3         1 


11-7T 


1 


(3-/2)    (/.- 1)^(1+,,^) 


/i-i 


-1] 


^[ 


(^-^^)  ^     (1  +  yo/ 


-li 
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For  n  =  3 


Pr{caplure\l)}  =  1 
-1 


1        -2A'i 


C2(3) 


z 

h=2 


KM 


-1 


^     (1  +  yo)' 


-  1])  +  ^o'E^"-! 


/!=3 


,4      10 


.M/2-1) 


(-•)'Tr7tr-"^^'7rT;5^-")^^^'^- 


.h{h-l){h-2)(^-^j^[ 

U  +  >0i  ^  /i=4 


(1+yo) 

6      7 


L^-.l-TT^l-^ 


(''-')'(i+yof 

/i!        /        1 


h=4 


3 

>'o 


+ 


C--^)      (l+y„) 


/!-2 


-1]-7T 


•[ 


1 


(1  +  yo) 


-1] 


1  -^Ar 


2^ 


(/2-4)!  V  (3-/^)    (i+y^f       (3-/2)  ■ 

"^  1  3         1 

= f 11  + -. 

13 


(3 -A)    (A-1)  '(!  +  ,„) 


■'(/-DItt-ttI 


.      -^ 11  +  ^ 


1 


/!      /     1         yo 


-i> 


U/-2)'(l+;,) 

1  A      -1 


-  1 


6      IS 


(l  +  vo) 

-^— I ^ 

('-2)      (!+>,)'■ 


7->l  +1^Z^"773 


+ 


(/-4)!   V(3-/)    (i  +  y^/       (3-/) 


-2 


Ij 


1]  + 


(3-/)      (/-I)    \\+y^t'  (^-^      ^ 


■[ 


(1  +  yo)' 


-11 
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For  n  =  4 


Pr{capture\4]  =  1  -  |  -^  e"^''°C2(4) 


z 

ft=2 


^.^(-^[  ^ 


^     (1  +  yo)' 


-1])  + 4X^.-1- 


Mh-\)(-n-\7[  ■■     ^  .._,  -i]  +  -r[       ^ 


4    il 
'0 


/!-l 


C-l)     (1  +  yo) 


.,„_,K.-2,(^,^ 


''    (1 +  >■»)' 


''  /l=4 


/i-2 


i]  +  7r^[        V.-.  -1] 


(^-1)    (i  +  vo) 


/!-l 


6     10 


/?! 


1 


3 

yo 


1  1  \  0    \^ 

"T'Trr;:^-";^!?^'^'-'  (,,_4):  v  {3-/0  (,+,„/ '  (3-y,)  ■ 


+ 


'  i^^^-ii    ' 


1 


e--^)        (l  +  y„) 


/I-2 


•[ 


1 


(l  +  ;'or 


--1] 


+  -T^  ^-'"'C3(4) 


(3 -A)    (/,-l)^(i^y^)^-i 


-ll  +  TT 


3         1 


(3-/z)    /i 

13 


(1  +  yo)' 


"'-'(T^'TTTtF-'^^'Titj-^-^S'^'- 


•'('-i)('-2)(-r:Vl        ' 


('-2)  '(I +  >■„)" 

6     )5 


U  +  TT^l     ,      '    ,-,   -1) 


(B.21) 


"-')'(i  +  yo) 


+  70)'        (^-'J 

1] I 2_   ^ 3        1 

/-2       'J      (3  _  A    r/-  n  ^  /,   ,      ^/-l       ^ 


(3-/)    (/-l)^(i+,^^) 


■[ 


1 


(1  +  yo) 


7-1] 


+  J_  ^-5^oV^(4) 


16 


(3-/)    / 

17 


Em  -h  1 77-^  - '  I  +  '^oE'^;-.  • 


y=4 


-^^     (1  +  yo/ 


;=5 


;/■(/■- 1)(  TT^  [ —^-TTT  -  1]  +  i  [ —^  -  1]    + -:r  E^y-2- 


(^■-1)    (i  +  yor 


•^'     (1  +  yo/ 


:/•(/■- 1)(/-- 2) 


-1  r  1  .,     .  2 


[ 


(^■-2)    (i  +  yo/ 


-r-il  +  7-^[ 


1 


-fiTTf;?-")**!'.- 


<7' 


(/■-i)^(i  +  yoy 

1 


-1] 


3 

yo        _^       3 


(/■-4)!  V  (3-7)   (l  +  vof       (3-7) 
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0--2)'(l  +  y„)'-^ 


-1| 


"J  ^  1 


ll  +  TT 


3         1 


(3-;)  ii-i)'  (I +  70/"'         ^^-^'^  >^'  ■ 


(l+Vo/ 


For  n  =  5 


Prlcapmr^lS)  =  1  -  |  -y  e-^^°C3(5) 


12 


1 


13 


EM-flTTT^-lll  +  ^o 


/=3 


(i+yo) 


E^'-i- 


/=4 


./(/-I) 


-1 


[ 


1 


./(/-!)(/- 2) 


(^-1)     (1  +  yo/ 
-1 


^-""^'TTw-'^-^S"'- 


('-2)     {l  +  y„) 


(1  +  Vo) 


4      14 


1 


-+'^-")*^5«.-. 


6      15 


^4i-i)     d  +  yo/-' 

/!        /       1 


1] 


rl      ^      3 


(/-4)!  V(3-/)    (i+,^y       (3-/) 
(/-2)      (1  +  yo)'-^         ^       (3-/)    (/-l)^(i  +  ,,/-i         ^       (3-/)    / 


(l+Vo) 


-  1 


2 


16 


17 


VAy-i[     ' 


ioZ'v-. 


4      1* 


*-  1)1  7^1-^-r-  n  +  il— ^-  11)  +4^ E^;-^- 


<^'-^)  ^  (i  +  yoy'       ^     -^  ^  (i  +  yo/ 


^■^■-^)^--K"fV^ 


1 


J-2 


1]  + 


^■-^)(i+yoy 


(i  +  yo/ 

1      r  1  ,A     ,       ^0      V,.  J-  /  1 


rr-1] 


-^    (i  +  yo/ 


-^lh^Z^-3  773^ 


"^°    r+     ^ 


^  '-'  (/•-4j!   V(3-;)    (l+yo/       (3-7)  • 


'       [-^^-l]-7T^7-^[-^^-l]+       ^        ' 


0--2)^(i  +  y„y-2 


(3-y)    (/■-!)  ^  (l+yo/"'  (^->^')   •^"  " 


(l  +  yo/ 


r-1] 


+  -y  e"^^°C5(5) 


20  21 

m=5        ^  u  +  yoJ  /         ^=6 
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■m{m  —  1) 


[ 


(^"-1)      (1  +  yor 


-1]  +  1[         ^ 


,4      22 


•m{m  —  l)(w  —  2) 


-1 


[ 


1 


(1  +  yo) 

2 


.  -1]    +    4 


M 


/^. 


w-2- 


m=i 


1 


(w  -  2)  '  (1  +  y  )"'-2        J  "^  ^m  _  n  ^  /I   ,  _   xw 


(^-1)     (1  +  yo) 


rr-1] 


1 


(1  +  yo) 


.6      23 
-^0     V 

36   l^' 


m\ 


1 


yoJ 

3 


+ 


/n  [  M   ,  -,  ^'^       ^U  +  1^  ^^""-^  (m  -  4)!  V  (3  -  /^)    (1  +  yo)'"    '    (3  -  w) 


1 


[ 


1 


(m-2)  ^  (l  +  yn)"-'       ^^ 


(3-m)    (m-  1)      (l+y^) 


vm-l 


-l]  +  7^ 


3         1 


(3-m)    ^  ■ 


•[ 


1 


(1  +  yo) 


m 


-1] 


And  finallv.  for  n  =  6 


Pr{capture\6]  =  I  — 


1      -aaI 


C,{6) 


12 


Y^^i' 


-1 


1 


/=3 


^      (l  +  y 


>'0i  '^  /=4 


■'''-'<  (^'7T7t?^-"^+'Trw- ")^^S"'- 


*-«'-<T/^'7T7tF 


''  1 

1]  +  77^^t[ 


I 


1 


^     (l  +  vo)' 


6     15 


-0-^Z-^- 


^^-^)     (1  +  yo)^-' 

/!       /       1 


-1] 


3 

yo 


+  Tr 


/=6 
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where  yo  is  the  preset  threshold  and  the  values  of  K^,  K,„  K„  Kj,  K„,  a.ndK^  are  given 
in  table  3. 
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